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Abstract 

Using events in which one of two neutral B mesons from the decay of an Y(AS) resonance is fully 
reconstructed, we set limits on the lifetime difference between the two neutral- -B mass eigenstates 
and on CP, T, and CPT violation in B°B° mixing. Both CP and non-CP eigenstates were obtained 
from the 88 million T(45) -> BB decays collected between 1999 and 2002 with the BABAR detector 
at the PEP-II asymmetric-energy P»-Factory at SLAC. We determine six independent parameters 
governing mixing (Am, AT/T), CPT /CP violation (Re-z, Imz), and CP /T violation (ImAcp, \q/p\), 
where Xcp characterizes B° and B decays to states of charmonium plus or K®. The preliminary 
results are 

sgn(ReAcp)Ar/r = -0.008±0.037(stat.) ± 0.018(syst.) [-0.084,0.068] , 
\q/p\ = 1.029±0.013(stat.) ±0.011(syst.) [ 1.001,1.057], 
(BbXcp/\Xcp\)Bjbz = 0.014±0.035(stat.) ± 0.034(syst.) [-0.072,0.101], 
Imz = 0.038±0.029(stat.) ± 0.025(syst.) [-0.028,0.104]. 

The values inside square brackets indicate the 90% confidence- level intervals. For lm\cp /\\cp\ an d 
Am we find values consistent with recent results from other analyses. These results are consistent 
with Standard Model expectations. 



Presented at the XXXVIII Rencontres de Moriond on 
Electroweak Interactions and Unified Theories, 
3/15-3/22/2003, Les Arcs, Savoie, France 

Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



The BABAR Collaboration, 

B. Aubert, R. Barate, D. Boutigny, J.-M. Gaillard, A. Hicheur, Y. Karyotakis, J. P. Lees, P. Robbe, 

V. Tisscrand, A. Zghiche 

Laboratoire de Physique des Particules, F-74941 Annecy-le-Vieux, France 

A. Palano, A. Pompili 
Universitd di Bari, Dipartimento di Fisica and INFN, 1-70126 Bari, Italy 

J. C. Chen, N. D. Qi, G. Rong, P. Wang, Y. S. Zhu 
Institute of High Energy Physics, Beijing 100039, China 

G. Eigen, I. Ofte, B. Stugu 
University of Bergen, Inst, of Physics, N-5007 Bergen, Norway 

G. S. Abrams, A. W. Borgland, A. B. Brcon, D. N. Brown, J. Button-Shafer, R. N. Cahn, E. Charles, 

C. T. Day, M. S. Gill, A. V. Gritsan, Y. Groysman, R. G. Jacobsen, R. W. Kadel, J. Kadyk, L. T. Kcrth, 
Yu. G. Kolomensky, J. F. Krai, G. Kukartscv, C. LeClcrc, M. E. Levi, G. Lynch, L. M. Mir, P. J. Oddone, 

T. J. Orimoto, M. Pripstcin, N. A. Roe, A. Romosan, M. T. Ronan, V. G. Shelkov, A. V. Telnov, 

W. A. Wenzel 

Lawrence Berkeley National Laboratory and University of California, Berkeley, CA 94720, USA 

T. J. Harrison, C. M. Hawkes, D. J. Knowles, R. C. Penny, A. T. Watson, N. K. Watson 
University of Birmingham, Birmingham, B15 2TT, United Kingdom 

T. Deppcrmann, K. Goctzcn, H. Koch, B. Lewandowski, M. Pelizaeus, K. Peters, H. Schmuecker, 

M. Stcinkc 

Ruhr Universitdt Bochum, Institut fur Experimentalphysik 1, D-44^80 Bochum, Germany 

N. R. Barlow, W. Bhimji, J. T. Boyd, N. Chevalier, W. N. Cottingham, C. Mackay, F. F. Wilson 
University of Bristol, Bristol BS8 1TL, United Kingdom 

C. Hearty, T. S. Mattison, J. A. McKenna, D. Thiessen 
University of British Columbia, Vancouver, BC, Canada V6T 1Z1 

P. Kyberd, A. K. McKemey 
Brunei University, Uxbridge, Middlesex UB8 3PH, United Kingdom 

V. E. Blinov, A. D. Bukin, V. B. Golubev, V. N. Ivanchenko, E. A. Kravchenko, A. P. Onuchin, 
S. I. Serednyakov, Yu. I. Skovpen, E. P. Solodov, A. N. Yushkov 

Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

D. Best, M. Chao, D. Kirkby, A. J. Lankford, M. Mandelkern, S. McMahon, R. K. Mommsen, W. Roethel, 

D. P. Stoker 

University of California at Irvine, Irvine, CA 92697, USA 
C. Buchanan 

University of California at Los Angeles, Los Angeles, CA 90024, USA 

H. K. Hadavand, E. J. Hill, D. B. MacFarlane, H. P. Paar, Sh. Rahatlou, U. Schwanke, V. Sharma 
University of California at San Diego, La Jolla, CA 92093, USA 



2 



J. W. Berryhill, C. Campagnari, B. Dahmes, N. Kuznetsova, S. L. Levy, O. Long, A. Lu, M. A. Mazur, 

J. D. Richman, W. Verkerke 

University of California at Santa Barbara, Santa Barbara, CA 93106, USA 

J. Beringer, A. M. Eisner, C. A. Heusch, W. S. Lockman, T. Schalk, R. E. Schmitz, B. A. Schumm, 
A. Seiden, M. Turn, W. Walkowiak, D. C. Williams, M. G. Wilson 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, CA 95064, USA 

J. Albert, E. Chen, M. P. Dorstcn, G. P. Dubois-Felsmann, A. Dvoretskii, D. G. Hitlin, I. Narsky, 

F. C. Porter, A. Ryd, A. Samuel, S. Yang 

California Institute of Technology, Pasadena, CA 91125, USA 

S. Jayatilleke, G. Mancinclli, B. T. Meadows, M. D. Sokoloff 

University of Cincinnati, Cincinnati, OH 4-5221, USA 

T. Barillari, F. Blanc, P. Bloom, P. J. Clark, W. T. Ford, U. Naucnberg, A. Olivas, P. Rankin, J. Roy, 

J. G. Smith, W. C. van Hoek, L. Zhang 

University of Colorado, Boulder, CO 80309, USA 

J. L. Harton, T. Hu, A. Soffcr, W. H. Toki, R. J. Wilson, J. Zhang 
Colorado State University, Fort Collins, CO 80523, USA 

D. Altenburg, T. Brandt, J. Brose, T. Colberg, M. Dickopp, R. S. Dubitzky, A. Haukc, H. M. Lacker, 
E. Maly, R. Muller-Pfefferkorn, R. Nogowski, S. Otto, K. R. Schubert, R. Schwicrz, B. Spaan, L. Wilden 

Technische Universitdt Dresden, Institut fur Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard, G. R. Bonneaud, F. Brochard, J. Cohen- Tanugi, Ch. Thicbaux, G. Vasileiadis, M. Verderi 
Ecole Polytechnique, LLR, F-91128 Palaiseau, France 

A. Khan, D. Lavin, F. Muheim, S. Playfer, J. E. Swain, J. Tinslay 
University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

C. Bozzi, L. Piemontese, A. Sarti 
Universita di Ferrara, Dipartimento di Fisica and INFN, 1-44100 Ferrara, Italy 

E. Trcadwell 

Florida A&M University, Tallahassee, FL 32307, USA 

F. Anulli, 1 R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, D. Falciai, G. Finocchiaro, P. Patteri, 

I. M. Pcruzzi, 1 M. Piccolo, A. Zallo 

Laboratori Nazionali di Frascati dell'INFN, 1-00044 Frascati, Italy 

A. Buzzo, R. Contri, G. Crosetti, M. Lo Vetere, M. Macri, M. R. Monge, S. Passaggio, F. C. Pastorc, 

C. Patrignani, E. Robutti, A. Santroni, S. Tosi 

Universita di Genova, Dipartimento di Fisica and INFN, 1-16146 Genova, Italy 

S. Bailey, M. Morii 
Harvard University, Cambridge, MA 02138, USA 
1 Also with Universita di Perugia, Perugia, Italy 



3 



G. J. Grenier, S.-J. Lee, U. Mallik 
University of Iowa, Iowa City, IA 52242, USA 

J. Cochran, H. B. Crawley, J. Lamsa, W. T. Meyer, S. Prcll, E. I. Rosenberg, J. Yi 
Iowa State University, Ames, IA 50011-3160, USA 

M. Davier, G. Grosdidicr, A. Hocker, S. Laplace, F. Le Diberder, V. Lepeltier, A. M. Lutz, T. C. Petersen, 

S. Plaszczynski, M. H. Schune, L. Tantot, G. Wormser 

Laboratoire de VAccelerateur Lineaire, F-91898 Orsay, France 

R. M. Bionta, V. Brigljcvic , C. H. Cheng, D. J. Lange, D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, CA 94550, USA 

A. J. Bevan, J. R. Fry, E. Gabathulcr, R. Garnet, M. Kay, D. J. Payne, R. J. Sloane, C. Touramanis 
University of Liverpool, Liverpool L69 3BX, United Kingdom 

M. L. Aspinwall, D. A. Bowerman, P. D. Dauncey, U. Egede, I. Eschrich, G. W. Morton, J. A. Nash, 

P. Sanders, G. P. Taylor 

University of London, Imperial College, London, SW7 2BW, United Kingdom 

J. J. Back, G. Bellodi, P. F. Harrison, H. W. Shorthouse, P. Strothcr, P. B. Vidal 
Queen Mary, University of London, El 4NS, United Kingdom 

G. Cowan, H. U. Flacchcr, S. George, M. G. Green, A. Kurup, C. E. Marker, T. R. McMahon, S. Ricciardi, 

F. Salvatore, G. Vaitsas, M. A. Winter 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, 

United Kingdom 

D. Brown, C. L. Davis 
University of Louisville, Louisville, KY 40292, USA 

J. Allison, R. J. Barlow, A. C. Forti, P. A. Hart, F. Jackson, G. D. Laffcrty, A. J. Lyon, J. H. Weatherall, 

J. C. Williams 

University of Manchester, Manchester M13 9PL, United Kingdom 

A. Farbin, A. Jawahery, D. Kovalskyi, C. K. Lae, V. Lillard, D. A. Roberts 
University of Maryland, College Park, MD 20742, USA 

G. Blaylock, C. Dallapiccola, K. T. Flood, S. S. Hertzbach, R. Kofler, V. B. Koptchev, T. B. Moore, 

H. Staengle, S. Willocq 

University of Massachusetts, Amherst, MA 01003, USA 

R. Cowan, G. Sciolla, F. Taylor, R. K. Yamamoto 
Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, MA 02139, USA 

D. J. J. Mangeol, M. Milek, P. M. Patcl 
McGill University, Montreal, QC, Canada H3A 2T8 

A. Lazzaro, F. Palombo 
Universitd di Milano, Dipartimento di Fisica and INFN, 1-20133 Milano, Italy 



4 



J. M. Bauer, L. Cremaldi, V. Eschenburg, R. Godang, R. Krocgcr, J. Rcidy, D. A. Sanders, D. J. Summers, 

H. W. Zhao 

University of Mississippi, University, MS 38611, USA 
C. Hast, P. Taras 

Universite de Montreal, Laboratoire Rene J. A. Levesque, Montreal, QC, Canada H3C 3J1 

H. Nicholson 

Mount Holyoke College, South Hadley, MA 01015, USA 

C. Cartaro, N. Cavallo, G. De Nardo, F. Fabozzi, 2 C. Gatto, L. Lista, P. Paolucci, D. Piccolo, C. Sciacca 
Universitd di Napoli Federico II, Dipartimento di Scienze Fisiche and INFN, 1-80126, Napoli, Italy 

M. A. Baak, G. Raven 

NIKHEF, National Institute for Nuclear Physics and High Energy Physics, 1009 DB Amsterdam, 

The Netherlands 

J. M. LoSecco 

University of Notre Dame, Notre Dame, IN 46556, USA 
T. A. Gabriel 

Oak Ridge National Laboratory, Oak Ridge, TN 31831, USA 

B. Brau, T. Pulliam 
Ohio State University, Columbus, OH 43210, USA 

J. Brau, R. Frcy, M. Iwasaki, C. T. Potter, N. B. Sincv, D. Strom, E. Torrence 
University of Oregon, Eugene, OR 91403, USA 

F. Colecchia, A. Dorigo, F. Galeazzi, M. Margoni, M. Morandin, M. Posocco, M. Rotondo, F. Simonetto, 

R. Stroili, G. Tiozzo, C. Voci 

Universitd di Padova, Dipartimento di Fisica and INFN, 1-35131 Padova, Italy 

M. Benayoun, H. Briand, J. Chauveau, P. David, Ch. de la Vaissiere, L. Del Buono, O. Hamon, 
Ph. Leruste, J. Ocariz, M. Pivk, L. Roos, J. Stark, S. T'Jampens 

Universites Paris VI et VII, Lab de Physique Nucleaire H. E., F-15252 Paris, France 

P. F. Manfredi, V. Re 
Universitd di Pavia, Dipartimento di Elettronica and INFN, 1-21100 Pavia, Italy 

L. Gladney, Q. H. Guo, J. Panetta 
University of Pennsylvania, Philadelphia, PA 19104, USA 

C. Angelini, G. Batignani, S. Bcttarini, M. Bondioli, F. Bucci, G. Calderini, M. Carpinclli, F. Forti, 
M. A. Giorgi, A. Lusiani, G. Marchiori, F. Martinez- Vidal, 3 M. Morganti, N. Neri, E. Paoloni, M. Rama, 

G. Rizzo, F. Sandrelli, J. Walsh 

Universitd di Pisa, Dipartimento di Fisica, Scuola Normale Superiore and INFN, 1-56121 Pisa, Italy 

2 Also with Universita della Basilicata, Potenza, Italy 

3 Also with IFIC, Instituto de Fisica Corpuscular, CSIC-Universidad de Valencia, Valencia, Spain 



5 



M. Haire, D. Judd, K. Paick, D. E. Wagoner 
Prairie View A&M University, Prairie View, TX 77446, USA 

N. Danielson, P. Elmer, C. Lu, V. Miftakov, J. Olson, A. J. S. Smith, E. W. Varncs 
Princeton University, Princeton, NJ 08544, USA 

F. Bellini, G. Cavoto, 4 D. del Re, R. Faccini, 5 F. Ferrarotto, F. Ferroni, M. Gaspero, E. Leonardi, 
M. A. Mazzoni, S. Morganti, M. Picrini, G. Piredda, F. Safai Tehrani, M. Serra, C. Voena 

Universitd di Roma La Sapienza, Dipartimento di Fisica and INFN, 1-00185 Roma, Italy 

S. Christ, G. Wagner, R. Waldi 
Universitat Rostock, D- 18051 Rostock, Germany 

T. Adyc, N. De Groot, B. Franek, N. I. Geddes, G. P. Gopal, E. O. Olaiya, S. M. Xella 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

R. Aleksan, S. Emery, A. Gaidot, S. F. Ganzhur, P.-F. Giraud, G. Hamel de Monchenault, W. Kozanecki, 
M. Langer, G. W. London, B. Mayer, G. Schott, G. Vasseur, Ch. Yeche, M. Zito 

DAPNIA, Commissariat a I'Energie Atomique/Saclay, F-91191 Gif-sur-Yvette, France 

M. V. Purohit, A. W. Weidemann, F. X. Yumiceva 
University of South Carolina, Columbia, SC 29208, USA 

D. Aston, R. Bartoldus, N. Berger, A. M. Boyarski, O. L. Buchmueller, M. R. Convcry, D. P. Coupal, 
D. Dong, J. Dorfan, D. Dujmic, W. Dunwoodie, R. C. Field, T. Glanzman, S. J. Gowdy, E. Grauges-Pous, 
T. Hadig, V. Halyo, T. Hryn'ova, W. R. Inncs, C. P. Jcssop, M. H. Kelsey, P. Kim, M. L. Kocian, 
U. Langenegger, D. W. G. S. Lcith, S. Luitz, V. Luth, H. L. Lynch, H. Marsiske, S. Menke, R. Messner, 
D. R. Mullcr, C. P. O'Grady, V. E. Ozcan, A. Perazzo, M. Perl, S. Pctrak, B. N. Ratcliff, S. H. Robertson, 
A. Roodman, A. A. Salnikov, R. H. Schindler, J. Schwiening, G. Simi, A. Snyder, A. Soha, J. Stelzer, 
D. Su, M. K. Sullivan, H. A. Tanaka, J. Va'vra, S. R. Wagner, M. Weaver, A. J. R. Weinstein, 
W. J. Wisniewski, D. H. Wright, C. C. Young 

Stanford Linear Accelerator Center, Stanford, CA 94309, USA 

P. R. Burchat, T. I. Meyer, C. Roat 
Stanford University, Stanford, CA 94305-4060, USA 

S. Ahmed, J. A. Ernst 
State Univ. of New York, Albany, NY 12222, USA 

W. Bugg, M. Krishnamurthy, S. M. Spanier 
University of Tennessee, Knoxville, TN 37996, USA 

R. Eckmann, H. Kim, J. L. Ritchie, R. F. Schwitters 
University of Texas at Austin, Austin, TX 78712, USA 

J. M. Izen, I. Kitayama, X. C. Lou, S. Ye 
University of Texas at Dallas, Richardson, TX 75083, USA 

4 Also with Princeton University, Princeton, NJ 08544, USA 

5 Also with University of California at San Diego, La Jolla, CA 92093, USA 



6 



F. Bianchi, M. Bona, F. Gallo, D. Gamba 
Universitd di Torino, Dipartimento di Fisica Sperimentale and INFN, 1-10125 Torino, Italy 

C. Borean, L. Bosisio, G. Delia Ricca, S. Dittongo, S. Grancagnolo, L. Lanceri, P. Poropat, 6 L. Vitale, 

G. Vuagnin 

Universitd di Trieste, Dipartimento di Fisica and INFN, 1-34127 Trieste, Italy 

R. S. Panvini 

Vanderbilt University, Nashville, TN 37235, USA 

Sw. Banerjee, C. M. Brown, D. Fortin, P. D. Jackson, R. Kowalewski, J. M. Roney 
University of Victoria, Victoria, BC, Canada V8W 3P6 

H. R. Band, S. Dasu, M. Datta, A. M. Eichenbaum, H. Hu, J. R. Johnson, R. Liu, F. Di Lodovico, 
A. K. Mohapatra, Y. Pan, R. Prcpost, S. J. Sckula, J. H. von Wimmersperg-Toeller, J. Wu, S. L. Wu, Z. Yu 

University of Wisconsin, Madison, WI 53706, USA 
H. Neal 

Yale University, New Haven, CT 06511, USA 



Deceased 



7 



1 Introduction and analysis overview 



The neutral B® meson system has two mass eigenstates with mass and total decay rate dif- 
ferences Am and Ar. While the mass difference has been measured recently with high precision 
[1, 2, 3, 4], only weak limits exist for the lifetime difference Atb = —AT/T 2 . Using the time- 
integrated mixing parameter Xd-, the CLEO Collaboration has set a limit of |Ar/r| < 80% [6]. A 
stronger constraint, | Ar/r| < 20% at 90% confidence-level, has been obtained by the DELPHI Col- 
laboration from a direct time-dependent study using flavor eigenstate events [7]. In the Standard 
Model, the ratio of the difference in the decay widths to the difference of the masses is proportional 
to ml/mf and thus quite small. Recent calculations of AT/T, including l/m^ contributions and 
part of the next-to-leading order QCD corrections within the Standard Model [5], find values of 
about —0.3%. The large data set available at the asymmetric B Factories provides the opportunity 
to reach much closer to the anticipated range for AT/T. 

The behavior of neutral B mesons is sensitive to CPT violation [8, 9]. The CPT theorem 
[10, 11], based on very general principles of relativistic quantum field theory, states that the triple 
product of the universal discrete symmetries C, P, and T represent an exact symmetry. The CPT 
symmetry remains to date the only combination of C, P, and T that is not known to be violated. 
However, the proof of the CPT theorem relies on locality, which could break down at very short 
distances. For instance, string theories are fundamentally non-local and therefore do not necessarily 
fulfill the conditions of the CPT theorem. Therefore it is possible, although perhaps unlikely, that 
CPT could break down. To date, the best tests have come from experiments in the neutral kaon 
system [12]. Bounds obtained so far in the B meson system [4, 13] are, however, mainly sensitive 
to the absorptive (lifetime) component of the Hamiltonian, where the small expected value of Ar 
suppresses the asymmetry effects. 

Violation of CP in the neutral B meson system may occur in mixing, in decay, or in the 
interference between mixing and decay. There is no fundamental way of assigning the source of CP 
violation observed in interference to either mixing or to decay. The standard phase choice puts CP 
violation in the mixing, but this is simply a convention. Other observable processes, however, can 
isolate CP violation due entirely to mixing. Similarly, mixing may intrinsically contain T violation 
or even CPT violation. It is these possibilities for the breaking of discrete symmetries in mixing 
itself that we address in this analysis. 

To measure the lifetime difference of the neutral .B-meson mass eigenstates and CP, T, or 
CPT violation we observe the time dependence of decays of neutral B mesons produced in pairs 
at the T(4S) resonance. The conventional mixing and CP analyses allow for exponential decay 
modulated by oscillatory terms with frequency Am. This neglects the difference between the decay 
rates Ar of the two mass eigenstates, which would introduce new exponential factors. CP, T, and 
CPT violation in the mixing of the neutral B mesons would modify the coefficients of the various 
terms involving exponential and oscillatory behavior. To detect these potential subtle changes 
requires precision measurements of the decays and thorough consideration of systematic issues. It 
also requires a more comprehensive treatment of the coherent decays of the mesons than has been 
conducted previously. 

The analysis is based on a total of about 88 million T(4S) — > BB decays collected between 1999 
and 2002 with the BABAR detector at the PEP-II asymmetric-energy B Factory at the Stanford 
Linear Accelerator Center. There, 9 GeV electrons and 3.1 GeV positrons annihilate to produce the 
BB pairs moving along the e~ beam direction (z-axis) with a Lorentz boost of (3j ps 0.56, allowing 
a measurement of the proper time difference At between the two B decays. In this analysis, one B 
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meson is fully reconstructed in a flavor (-Bfl av ) or CP (Bqp) eigenstate (generally denoted as B rec ). 
The remaining charged particles in the event, which originate from the other B meson (B taig ), are 
used to identify its flavor as B° or B°. The time difference At = t rec — t tag ~ Az/(3^/c is determined 
from the separation Az of the decay vertices for the fully reconstructed B candidate and the tagging 
B along the boost direction. 

A single maximum-likelihood fit to the time distributions of tagged and untagged, flavor and 
CP eigenstates determines six independent parameters (see Sec. 2) governing mixing (Am, AT/T), 
CPT/CP violation (Rez, Imz) and CP /T violation (ImAcp, \q/p\), where Xcp is the traditional 
variable characterizing the decays of neutral B mesons into final states of charmonium and a K® or 
K® . The parameters ImAcp and Am are used only as a cross-check with the BABAR sin 2/3 analysis 
[14] and previous Am results [1, 2, 3, 4]. 

The analysis has several challenges. First, the tagged B and the fully reconstructed B decays are 
correlated and interference between allowed and doubly-CKM-suppressed (DCKM) decays cannot 
be neglected. Second, tagging incorrectly assigns the flavor with a certain mistag probability. Third, 
the resolution for At is comparable to the B lifetime and asymmetric for positive and negative At. 
This asymmetry must be well understood lest it be mistaken for a fundamental asymmetry we seek 
to measure. Fourth, possible direct CP violation in the Bqp sample can be a competing source 
of fake effects and must be parameterized appropriately. Finally, we have to account possible 
asymmetries induced by the differing response of the detector to positive and negative particles. In 
resolving these issues we rely mainly on data. 

The paper is organized as follows. In Sec. 2 we present a general formulation of the time- 
dependent decay rates of B°B° pairs produced at the T(4S") resonance, including effects from 
the lifetime difference, possible CPT violation, and interference effects induced by doubly-CKM- 
suppressed decays. We derive the expressions for B decays to final states with flavor and CP 
eigenstates. In Sec. 3 we describe the BABAR detector. After discussing the data sample in Sec. 4, 
we describe the 6-quark tagging algorithm in Sec. 5. Sec. 6 is devoted to the description of the 
measurement of Az and the determination of At and its resolution function. In Sec. 7 we describe 
the unbinned log- likelihood function and the assumptions made in the nominal fit. The results of 
the fit are given in Sec. 8. Cross-checks are discussed in Sec. 9 and systematic uncertainties are 
summarized in Sec. 10. The results of the analysis are summarized and discussed in Sec. 11. 



2 General time-dependent decay rates from T(45) — > B°B° 

The neutral B meson system can be described by the effective Hamiltonian M — iT/2, where 
M and F are two-by-two hermitian matrices describing, respectively, the mass (dispersive) and 
lifetime (absorptive) components. If either CP or CPT is a good symmetry, then Mn = M22 
and Tn = T22, with the index 1 indicating B° and 2 indicating B°. If either CP or T is a good 
symmetry, T12/M12 is real. This condition does not depend on the phase convention chosen for the 
B° and B°. The masses \i± and decay rates j± of the two eigenstates are 



"± = M± " \l± = M - l -T ± J (m 12 - l -T 12 ^j (V* 2 - ^r* 2 ) + (dM - l -5T^ , (1) 



where we define 



M11 + m 22 Tn + r 2 2 M11 - m 22 Tn - r 2 2 

M EE , r EE , 5M EE — , <5r EE . (2) 
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Neglecting CPT violation, and anticipating that the lifetime difference is small compared to the 
mass difference, we have 



Am = 2IM12I 



Ar = 2|Mi 2 |Re(r 12 /Mi 2 ) 



(3) 



Here we have taken Am to be the mass of the heavier state minus the mass of the lighter. Thus 
Ar is the decay rate of the heavier state minus the decay rate of the lighter and its sign is not 
known a priori. 

The light and heavy mass eigenstates of the neutral S-meson system may be written 



\Bl) 
\Bh) 



p\B°)+q\B } 
p\B°) - q\B°) 



(4) 



where 



q 
p 



The magnitude of q/p is very nearly unity: 



I M 



12 



2 1 12 



Mi 2 - i r 12 



(5) 



1 — Im- 



12 

"M12" 



(6) 



In the Standard Model, the CP- and T- violating quantity \q/p\ 2 — 1 is small not just because |ri 2 | is 
small, but additionally because the CP-violating quantity Im(ri 2 /Mi2) would vanish if the u- and 
c-quark mass were the same. CP violation is not possible in mixing if two of the quark masses (for 
quarks of identical charge) are identical because we could redefine them so one quark did not mix 
with the other two. The remaining two generations would be inadequate to support CP violation. 
The result is that Im(ri2/Mi 2 ) is suppressed by an additional factor m\lm\ ~ 0.1. When the 
remaining factors are included, the result is |Im(ri2/Mi 2 )| < 10 -3 . 

CPT violation in mixing can be described conveniently by the phase-convention independent 
quantity 



z = 



5M - § <5r 



SM - § 5T 



sj (Mi 2 - § T 12 ) (M* 2 - § T* 12 ) + (6M-i ST) 2 \ ( Am " § Ar ) 



(7) 



States that begin as purely B° or B° will oscillate and after a time t will be mixtures 



|5° phys (t)> = (g + (t) + z ■ g-(t)) \B°) - VT 



z 2 --9-(t)\B°) 



I -D phys 



(*)) = (g + (t) ~ z ■ g-(t)) \B°) - v / T^- - q9 -{t) \B°) , 



(8) 



where we have introduced 



3 ±{t) = \{e- i "+ t ±e- i "- t ) 



(9) 



At the T(4S) resonance, neutral B mesons are produced in coherent pairs. If we subsequently 
observe a final state f± at time to = and another state / 2 at some other time t, either positive 
or negative, we cannot in general know whether f\ came from the decay of a B° or a B° and 
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similarly for the state / 2 . If ^1,2 an d A 12 are the amplitudes for the decay of B° and B° to the 
states f\ and f 2 , then the overall amplitude when t > is given by 

A = a + g + {t) + a-g-{t) , (10) 

where 

a+ = -A 1 A 2 + A X A 2 

+ z[A{A 2 +^lA 2 ] . (11) 



a_ = \A — z 2 



V -A X A 2 - ^-A X A 2 



P 



Using the relations 



\g±(t)\ 2 = - e~ t/TB [cosh(AIY/2) ± cos(Amt)] (12) 



dN _ r|; 
— — oc e 
dt 



and 

g* + {t) g-(t) = -\z~ ttTB [sinh(Art/2) + i sin(Amt)] , (13) 
with tb = 1/r, we find the decay rate, which in fact is correct for t positive or negative, 

!{-c + cosh(AR/2) + ^c_cos(Amt) - Res sinh(AR/2) +Ims sin(Amt)} (14) 
where 

c ± = l a +l 2 =t l a -l 2 ; s = o+o- • (15) 

Now let us take f\ = ftag to be the state that is incompletely reconstructed and which provides 
the tagging decay, and / 2 = free the fully reconstructed state (flavor or CP eigenstate). Because 
the tagging algorithm is imperfect, we may incorrectly identify the flavor of the decaying B meson. 
This can be accounted for by incoherently combining correct and incorrect tags. A more subtle 
problem arises because there may be a basic ambiguity: the state ftag may result from interference 
between decay from a B° and decay from a B°. We consider first the simpler situation where there 
is no underlying ambiguity. 

If the tag is a B°, we display this explicitly writing Aj^ / 0, A-^ = 0. We define 

q A rec 

*rec — -—. , UDJ 

P s±rec 

which is independent of phase conventions for the B° and B° states. Dropping terms of order z 2 , 
we find a decay rate 

^L(tag = B°) oc \A^\ 2 \A rec \ 2 e'^ j \ [l + | A rec - z\ 2 ] cosh(An/2) 

+^ [1 - \Kec - z\ 2 ] cos(Ami) 

-Re(-A rec + z) sinh(An/2) + Im(-A rec + z) sin(Amt) \ . (17) 
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Correspondingly, if f tag is tagged as a B°, A tag = 0, A tag / 0, and we have 
[N (tag = B°) oc \A tag \ 2 \A rec \ 2 



dt 



P 



[|A rec | 2 + |l + zA rec | 2 ] cosh(Art/2) 



+\ [\Kec\ 2 - |1 + z\ rec \ 2 ] cos(Amt) 

+Re[A; ec (l + Xrecz)} sinh(Art/2) - Im[A; ec (l + X rec z)} sin(Ami) j . 

(18) 

The normalizations are identical in Eqs. (17) and (18). 

We consider several scenarios. If free is & CP eigenstate, then |A(7p| — Pcp\q/p\i with Rqp — 
\Acp/Acp\- If all the weak decay mechanisms have the same weak phase, Rcp = 1- This is expected 
for final states like J/ipK®, where indeed measurements show |Acp| ~ 1 [14]. Dropping quadratic 
terms in z and Ar we have 



dN 
~dt 



(tag = B°;rec = CP) oc \A^\ 2 \ A CP \ 2 e- m j 1 [l + \\ CP \ 2 - 2ReAcpRez - 2ImAcpIi 

+- [1 - \\cp\ 2 + 2ReA C pRez + 2ImA C pImz] cos(Amt) 



2 

+ReAcp sinh(Art/2) - Im [X C p - z] sin(Amt) \ , (19) 



H /V n 2 

— (tag = B°;rec = CP) oc \A tag \ 2 \A CP \ 2 V - 



e ~ m |^ t 1 + \ X CP? + 2ReA C pRez - 2ImAcpImz] 



+ \ [|Acp| 2 - 1 - 2ReA C pRez + 2ImAcpImz] cos (Ami) 
+ReA C p sinh(Art/2) + Im [\ C p - \Xcp\ 2 z] sin(Ami) | . (20) 

Data from directly related final states like J/ipK®, with tjcp = — 1, and JfyK®, with rjcp = +1, 
where r\cp is the CP eigenvalue of the final state, can be combined by assuming that they are 
identical, except for an overall sign in Acp- 

We assume that the decays of flavor eigenstates are dominated by a single weak mechanism, so 
that |^4flav| = l^fi^l) l^flavl = I^Havl- This will enable us to relate the four possibilities that arise 
from the tag and reconstructed state being either B° or B°. When the fully reconstructed meson 
f rec is a flavor eigenstate, |A rec | is either very small or very large. If it appears to come from a B°, 
then |Afl av | <C 1 and we have 
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^(tag = B°; rec = B°) oc | AaglVflav| 2 e~ r|t| \\ + \ cos(Amt) + Im(-A flav + z) sin(Amt) \ , 



dN 



(tag = B°;rec = B°) oc \A tag \ 2 \A aav 



P 



(21) 



e r '*' <--- cos(Amt) + ImAfl av sin(Amt) > . 



(22) 



Conversely, if the fully reconstructed state is nominally a B°, | Ag^| 3> 1 and 



dN 



— «<1 = B°;rec = B°) oc |^tag| 2 |^| 2 |A BS7 | 2 e- r l^ ^ - \ cos( A,/;/) + L,nA (lav sin(A//,/ j 



OC lAaglVflavl 2 



-r|*|. 



1 1 

2 ~ 2 



cos (Ami) + ImAflav sin(Amt) > , 



(23) 



^r(tag = B°;rec = B°) oc |^ta g | 2 |^| 2 



I A: 



+ - cos(Ami) 



— lm[z + Aflav] sin(Ami) 



oc |Aag| 2 |^4flav| 2 e r| *' \\ + \ cos(Amt) - lm[z + A flav ] sin(Amt) 



(24) 



where Afl av = 1/Ag^. We see that the overall normalization in the mixed final states has a factor 
of \q/p\ 2 or \p/q\ 2 relative to the unmixed final states. 

Doubly-CKM-suppressed decays, such as B° — > D + ir~, occur at a rate roughly \V* b V c d/V* b V u d\ 2 ~ 
(0.02) 2 , and can be ignored. However, interference between favored and suppressed decays are 
reduced by a factor of approximately 0.02. For decays reconstructed in a final state as appar- 
ent B° mesons, we anticipate |Afl av | « 0.02|g/p|, while for reconstructed B° mesons, |Afl av | ~ 
0.02|p/g|. In principle, every hadronic final state has a different Afl av that can be written as 



Afla 



I Aflav I e" 



where 4> s is the strong phase and 4> w is the weak phase. Provided that a 



single mechanism contributes to the allowed and suppressed decays, Afl av and Afl av have the same 
strong phase but opposite sign weak phase, and the magnitudes are the same up to a relative \p/q\ 2 
factor. For simplicity, in our analysis we take the DCKM effect to be equal for all reconstructed 
flavor states. 

Eqs. (17)-(24) show that while ImAcp, Imz, \q/p\ and Rqp are unambiguously determined, Rez 
appears only in the product ReAcpRez. Similarly, Ar cannot be determined separately from ReAcp 
since there is an ambiguity in ReAcp: ReAcp = ±\/|Acp| 2 — (ImAcp) 2 - As a result, the parameters 
which can actually be determined by the analysis are sgn(ReAcp)Ar/r, \q/p\, (ReAcp/|Acp|) Rez, 
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Table 1: Dominant sensitivity of physical parameters to the distributions measured with the fully 
reconstructed flavor and CP states. The flavor sample is much larger than the CP sample. 



-Bflav B C p 
Variable t-even i-odd t-even t-odd 

sgn(ReAcp)Ar/r x 
\q/p\ x x 

(ReA C p/|A C p|)Rez x 
Imz x 

lm\ C p/\Xcp\ x 
R C p x 
Am x 



Imz, lm\cp/\\cp\, Rcp, Am and Tp = 1/T. Both CP eigenstates and flavor eigenstates are needed 
for the analysis, as shown in Table 1. The sensitivity to (ReAcp/|Acp|) Rez and ImAcp/|Acp| is 
provided by the CP-eigenstate events P>cp, for which the t dependence is even for the former and 
odd for the latter. The sample contributes marginally because it lacks explicit dependence on 
ImAcp/|Acp| and the dependence with Rez is scaled by the sinh (Art/2) term, which is small for 
small Ar. In contrast, \q/p\ and Imz (and Am) are completely dominated by the large statistics 
-Bflav sample, for which the t dependence is even for the former and odd for the latter. For small 
values of AT/T, the determination of AT/T is dominated by the B>cp sample, in spite of the 
relatively small statistics compared to the f?fl av sample. This is due to the even cosh(Art/2) 
dependence (Ar 2 to first order) of the flavor sample, while the CP sample has a non-vanishing 
odd sinh(Art/2) (Ar to first order) dependence. The contribution of sinh(Art/2) is the same 
for both B° and B° tags, so untagged events may be included as well. The Bcp sample is also 
sensitive to the sign of AT/T (up to the sign ambiguity from ReAcp). Overall, the combined use 
of the -Bflav an d Bcp samples provides maximal sensitivity to the physical parameters, since they 
are determined either from different samples, or from different t dependencies. Small correlations 
are induced by the detector resolution. 

Doubly-CKM-suppressed decays occur on the tagging side as well, if the decay is non-leptonic. 
Consider, for example, the contribution of doubly-CKM-suppressed decays when the tagging decay 
is ostensibly a B°. To first order in A tag = (q/p)Ata,g/A ta g, z and AT/T, the new contribution to 
the decay rate is 

dN 



■-(tag = B°; DCKM contribution) oc | ^tag 1 2 1 ^rec 1 2 



2 

„-r|t| 



(25) 



x \ ^ [-4ReA tag ReA rec ] + ^ [-4ImA ta gImA rec ] cos(Ami) + ImA ta g(|A rec | 2 - 1) sin(Amt) 



This shows that if the reconstructed state is a flavor eigenstate, the DCKM effect in tagging is 
negligible except in the sin(Amt) term. Conversely, if the reconstructed state is a CP eigenstate 
with |A rec | ~ 1, the DCKM effect is confined to the terms even in t. Because terms quadratic in 
Atag can be ignored, the combined effect of DCKM on tagging can be incorporated in one value of 
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A tag and one of A tag = 1/A^. 

3 The BABAR detector 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The BABAR detector is described in detail elsewhere [16], so here we provide only a brief 
description of the apparatus. 

Surrounding the beam-pipe is a five-layer silicon vertex tracker (SVT), which provides precise 
measurements of points along the trajectories of charged particles as they leave the interaction re- 
gion. This allows track reconstruction, even for some particles with momentum less than 120 MeV/ c. 
Outside of the SVT is a 40-layer drift chamber (DCH) filled with an 80:20 helium-isobutane gas 
mixture, chosen to minimize multiple scattering. The DCH measurements provide charged-particle 
tracking and determination of momenta through track curvature in the 1.5-T magnetic field gener- 
ated by the superconducting coil. The DCH also provides dE/dx energy- loss measurements, which 
contribute to charged-particle identification. Surrounding the drift chamber is a novel detector of 
internally reflected Cerenkov radiation (DIRC), which provides charged-particle identification in 
the barrel region. Outside of the DIRC is a CsI(Tl) highly segmented electromagnetic calorimeter 
(EMC), which is used to measure the energy of photons, to provide electron identification, and to 
detect neutral hadrons through shower shapes. Finally, the flux return of the superconducting coil 
surrounding the EMC is instrumented with resistive plate chambers interspersed with iron (IFR) 
for the identification of muons and neutral hadrons. 

A detailed Monte Carlo program based on the GEANT4 [17] software package is used to simulate 
the BABAR detector response and performance. The agreement between data and simulation is very 
good [16]. 

4 Data samples and B meson reconstruction 

We have selected events where one of the B mesons is completely reconstructed in either a 
flavor (i?fl av ) or CP (Bcp) eigenstate, using the same criteria used for the BABAR hadronic Am 
[1, 18] and sin 2/3 measurements [14]. The decay modes used for the flavor sample, the CP sample, 
and a control sample are displayed in Table 2. 

We select Sfl av and Bqp candidates by requiring that the difference AE between their energy 
and the beam energy in the center-of-mass frame be less than three standard deviations from zero. 
The AE resolution ranges between 10 and 50 MeV depending on the decay mode. For -Bfl av and Bqp 
modes involving K® (B CPK o), the beam-energy substituted mass must be greater than 5.2GeV/c 2 . 

The beam-energy substituted mass is given by toes = \J (\s + p i ■ ps) 2 /Ef — p 2 B , where s is the 
square of the center-of-mass energy, Ei and p» are the total energy and the three-momentum of the 
initial state in the laboratory frame, and ps is the three- momentum of the B candidate in the same 
frame. In the case of decays to JfyK®, the K\ direction is measured but its momentum is only 
inferred by constraining the mass of the J/ipK® candidate to the known B° mass. As a consequence 
there is only one parameter left to define the signal region, which is taken to be \AE\ < 10 MeV. 

The purities are determined from fitting the data to the tties (-Bflav and B CPK o modes) or AE 
(B CPK o mode) distributions [18]. Figure 1 shows the wes distribution for the -Bfl av and B CPK o 
samples and the AE distribution for the B CPK o candidates, before the vertexing requirements 
(see Sec. 6). The combinatorial background in the uies distributions is described by an empirical 
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Table 2: The flavor, CP, and control samples used in this analysis. Charged and neutral flavor 
eigenstate decay modes imply also their charge conjugate. 

Samples Decay modes 



-Sflav 

B° D*-ir+(p+,a + ) 



D*- -► D°tt- 



5° -► J D-^+(p+,a+) 
5° -► J/i>K*° 



















J/0 — > e + e , 






^^7T + 7T-,7r 7r 












^(25) -> e + e~,fi + fj,~,J/ip vr+vr" 


















Xci -> J/0 7 






K° s -► 7T+^- 




-► 





Control 




5+ - 7J W V 


77° - 73%° 


5+ -► j/v^k + 








5+ - XdK+ 




B+ -► J/0;K*+ 





phase-space model [18] and the signal with a Gaussian distribution. The combinatorial background 
consists of continuum and BB sources, and has a time structure with both prompt and non-prompt 
components. A small correlated background due to other B decays (not shown) also peaks at the 
B mass. The background in the J/ipK® channel receives contributions from other B decays with 
real J/0 mesons in the final state, and combinatorial sources. 

After completely reconstructing one B meson, the rest of the event is analyzed to identify the 
flavor of the opposite B and to reconstruct its decay point, as described in Sees. 5 and 6. 

Using exactly the same requirements, we analyze GEANT4-simulated samples of generic BB 
and signal events to check for any biases in the procedure or extracted parameters. The Monte 
Carlo samples are also used to assess detector response and to estimate some background sources. 
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Figure 1: Distributions for -Bfl av and Bcp candidates before vertexing requirements: a) tues fo r 
-Bflav states; b) toes fo r B° — > J/ipK®, ip(2S)K®, XciK® final states; and c) Ai? for the final state 
5° -> J/0-fifg. 

The values of the physics parameters assumed in the simulations are similar to those measured in 
the data. We used additional samples with significantly different values to check the reliability of 
the analysis in other regions of the parameter space. 

5 Flavor tagging 

The tracks that are not part of the fully reconstructed B meson are used to determine whether 
the -Btag was a B° or B° when it decayed. This determination cannot be done perfectly. If 
the probability of an incorrect assignment is w, an asymmetry that depends on the difference 
between B° and B° tags will be reduced by a factor D = 1 — 2w, frequently called the dilution. A 
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neural network combining the outputs of physics-based algorithms is used to take into account the 
correlations between the different sources of flavor information and to assign the event to one of 
five mutually exclusive tagging categories. The dilution for each category is determined from data. 
Grouping tags into several categories, each with a relatively narrow range in mistag probability, 
increases the overall power of the tagging. 

We group together events of similar character to make it possible to study systematic effects. 
Events with an identified primary electron or muon and a supporting kaon, if present, are assigned 
to the Lepton category. The KaonI category contains events with an identified kaon and a soft- 
pion candidate with opposite charges and similar flight direction. Soft pions from D* + decays are 
selected on the basis of their momentum and direction with respect to the thrust axis of B tag . 
Events with only an identified kaon are assigned to the KaonI or KaonI I category depending on the 
estimated mistag probability. Events with only a soft-pion candidate are assigned to the KaonI I 
category as well. The remaining events are assigned to the Inclusive or UnTagged category based 
on the estimated mistag probability. The UnTagged tagging category has a mistag rate near 50%, 
and therefore does not provide tagging information but it increases the sensitivity to the lifetime 
difference and allows the determination from the data of the detector charge asymmetries, as 
described in Sec. 7. The tagging efficiencies t a , defined after vertexing cuts, for the five categories 
are measured from the data and summarized in Table 4 of Sec. 8. This tagging algorithm is identical 
to that used in Ref. [14]. 

The mistag probabilities appear separately for B° and B° tags in each tagging category. If we 
define w a = (wg + w^ )/2 and Aw a = w B0 — w^ , we have that 

w b°b° =w a ± A ^72 • (26) 

A correlation between the mistag rate and the At uncertainty estimated event-by-event (discussed 
in Sec. 6) is observed in the Monte Carlo simulation for kaon based tags [2, 18]. For a At uncertainty 
less than 1.4 ps, this correlation is found to be approximately linear: 

w a = wg + w2 ope CTAt. (27) 
All signal mistag parameters, Wq, w® lope and Aw a , are floated in the global fit (11 in total since 
u>siopc° n is assumed to be zero), and their results can be found in Table 4 in Sec. 8. 

6 Decay time measurement and At resolution function 

The time interval At = t rec — t tag between the two B decays is calculated from the measured 
separation Az between the decay vertex of the reconstructed B Tec meson and the vertex of the -B tag 
meson along the z-axis, using the known boost of the T(4<S) resonance in the laboratory, (3~f = 0.56. 
The method is the same as described in Sec. V in Rcf. [18]. 

An estimated error a^t on At is calculated for each event. This error accounts for uncertain- 
ties in the track parameters from the SVT and DCH hit resolution and multiple scattering, our 
knowledge of the beam spot size and most of the effects from the average B flight length in the 
vertical direction as well as the Az to At conversion. However, it does not account for errors due to 
mistakes of the pattern recognition system, wrong associations of tracks to vertices, misalignment 
within and between the tracking devices, inaccuracies in the modeling of the amount of material in 
the tracking detectors and in our knowledge of the beam spot position and size, and the absolute 
z scale. We use parameters in the At resolution function, extracted from the data in the global 
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fit (Table 5 in Sec. 8), to absorb most of these effects. Remaining systematic uncertainties are 
discussed in detail in Sec. 10. 

We use only those events in which the vertices of the B rec and B tag are successfully reconstructed 
and for which \At\ < 20 ps and the At uncertainty (<JAt) is less than 1.4 ps, as used in Ref. [1]. The 
fraction of events in data satisfying these requirements is about 85%. From Monte Carlo simulation 
we find that the reconstruction efficiency does not depend on the true value of At. Excluding 0.3% 
of the events that are poorly measured, the rms vertex resolution is about 160 //m (1.0 ps). Because 
this resolution is poorer than that for the completely reconstructed vertex, the overall Az resolution 
is dominated by the resolution of the B tag vertex. 

To model the At resolution we use the sum of three Gaussian distributions (called core, tail 
and outlier components) with different means and widths: 



lZ(5t,aAt; tf) = (1 - /tail - /out)^G(<ft; S cole aAt, S corc aAt) + 
/tail^G^; ^tail^At, S tai \a At) + 

/out^G(<ft;<5out)<7but) (28) 

with 

h G (5t;5,a) = -;L- exp(-(<5t - 8) 2 /(2o 2 )) , (29) 

where q represents the collection of parameters needed to describe the Gaussians and St = At — At true 
represents the reconstruction error. The vertex reconstruction provides an event-by-event estimate 
for \St\, namely a At- We incorporate the last Gaussian in Eq. (28) without reference to a At since 
the outlier component is not expected to be well described by the estimated error. The first two 
Gaussian components allow two independent scale factors, Score and Staib to accommodate an over- 
all underestimate (> 1) or overestimate (< 1) of the estimated errors. The core and tail Gaussian 
distributions are allowed to have non-zero means (8 core and 5 ta ih respectively) to account for residual 
biases due to daughters of long-lived charm particles included in the B tag vertex. Separate means 
are used for the core distribution of each tagging category. These means are scaled by a At to ac- 
count for a correlation observed in Monte Carlo simulation between the mean of the St distribution 
and a At [2, 18]. This correlation is found to be approximately linear for a a At less than 1.4 ps. The 
non-zero means of the resolution functions introduce an asymmetry into the otherwise symmetric 
At distributions. The outlier Gaussian has a global width (cr ou t) an d offset (<5 ou t), and it accounts 
for less than 0.3% of the reconstructed vertices. The parameters /tail and / ou t in Eq. (28) represent, 
respectively, the fractions of the tail and the outlier Gaussian components. In simulated events, we 
find no significant differences between the At resolution function of the -Bfl aV) B CPK o and B CPK o 
samples. This is expected, since the B ta , g vertex precision dominates the At resolution. Hence, the 
same resolution function is used for all modes. Residual differences are taken into account in the 
evaluation of systematic errors, as described in Sec. 10. 

We find that the three parameters describing the outlier Gaussian component are largely cor- 
related among themselves and with other resolution function parameters. Therefore, we fix the 
outlier bias <5 ou t and width <7 ou t to and 8 ps, respectively, and vary them in a wide range to 
evaluate systematic uncertainties. The resulting signal resolution function is described by a total 
of 12 parameters, 

-» / n cLepton cKaonl rKaonll ^Inclusive rUnTagged f X Q f X rr \ 

y — \ '-'core; °core > °core > °core > °corc > °corc j /tail) °tail) '-'tail; /out) "out; Oout } > 
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10 of which are floated in the final fit. Their final values can be found in Table 5 in Sec. 8. 

As a cross-check, we use an alternative resolution function that is the sum of a single Gaus- 
sian distribution (centered at zero), the same Gaussian convolved with a one-sided exponential to 
describe the core and tail parts of the resolution function, and a single Gaussian distribution to 
describe the outlier component [2]. The exponential component is used to accommodate the bias 
due to tracks from charm decays on the Stag side. The exponential constant is scaled by a A t to 
account for the previously described correlation between the mean of the 5t distribution and a A f In 
this case, each tagging category has a different core component fraction and exponential constant. 



7 Likelihood fit method 

We perform a single unbinned maximum-likelihood fit to all i?fl av , B CPK o and B CPK o events. 
Each event is characterized by its assigned tag category, a £ {Lepton, KaonI, Kaonll, Inclusive, 
UnTagged}; its tag-flavor type, tag = B°,B° (unless it is untagged); its reconstructed event type, 
rec = B°, B°, CPK®, CPK®; the values of At and aAt', and a variable £, either tties or AE, used 
to assign the probabilities that the event is signal or background. The underlying distributions 
depend on whether the event is signal or any of a variety of backgrounds (together specified by j), 
on the tag category (a), on the tag flavor (tag), and on the reconstructed final state type (rec). 
The contribution of a single event to the log-likelihood is 



In 



E 



^i(onti rec (At,a At ) 



(30) 



For a given reconstructed event type rec and tagging category a, J>ec(C) gives the probabilities 
that the event belongs to the signal or any of the various backgrounds indicated by j. Each 
such combination has its own PDF H^ grec (At, aAt), which depends as well on the particular tag 
flavor, tag. This distribution is a convolution of a tagging-category-dependent time distribution, 
H^ rec (At tr uc), with a At resolution function 



/+oo 
d(At truc )^(At - At tIUC ,a At ;q^)H^ rec (At tTUC ) (31) 
-oo 

where 

^ta<7,rec(Aitrue) = r rec \^tag 0- ~ w tag)Hag,rec(^truc) 

Here h J tag rec (Att rue ) represents the time dependence d-ZV/dAttme described in Sec. 2, with At tr uc = t, 
using full expressions rather than expansions in Ar and z. We indicate by u>^^— the mistag frac- 
tions for category a and component j. The index tag denotes the opposite flavor to that given by 
tag. For events falling into tagging category UnTagged, w ^g[tag = efficiency t"^ g is the 

probability that an event whose signal/background nature is j and whose true tag flavor is tag will 
be assigned to category a, regardless of whether the flavor assigned is correct or not. The efficiency 
riec is the probability that an event whose signal/background nature is indicated by j and whose 
true character is rec will, in fact, be reconstructed. 
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7.1 PDF normalization 



Every reconstructed event, whether signal or background occurs at some time At true , so 



j — < 



d(At true )/i^ s rec (Attrue) = 1, for all rec, tag and j . (33) 
Moreover, every event is assigned to some tag category (possibly UnTagged), thus 

£C^ = 1- (34) 

a 

It follows then that the normalization of H££ rec (At true ) is 

/+0O 
d(Ai true )<£ rec (Ai true ) = ri ec . (35) 
-oo 



EE 



a tag 

In this analysis the nominal normalization of Ti.^ grec {At, cta*) is the same as (Att rU e) 
(asymptotic normalization), but fits with normalization in the interval [—20,20] ps were also per- 
formed as a cross-check to evaluate possible systematic effects. 

7.2 Efficiency asymmetries 

For each signal or background, j, there is an average reconstruction efficiency, r 7 = (r^„+r^ )/2. 
These average efficiencies are ultimately absorbed when we define fractions of reconstructed events 
falling into the different signal and background classes. In contrast, because all events fall into 
some tagging category (including UnTagged), the averages t a >i = (t^ + t^)/2 are meaningful, 
and for j=signal, and a = UnTagged the result plays an important role. The asymmetries in the 
efficiencies, 



v 3 = 



+ot,3 _ ,a,j 

^ = B °. S °. , (36) 

Z B° Z B° 

need to be determined because they might otherwise mimic fundamental asymmetries we seek to 
measure. In Appendix A we illustrate how the use of the untagged sample makes it possible to 
determine the asymmetries in the efficiencies. Note that asymmetries due to differences in the 
magnitudes of the decay amplitudes, lAflavl / l^iwl an£ f l-^tagl l-^tagl> cannot be distinguished 
from the asymmetries in the efficiencies, thus are absorbed in the v, /j, a parameters. 

We fix the average tagging efficiencies t a ' J to estimates determined by counting the number 
of events falling into the different tagging categories, for each decay channel separately and in all 
the C range, after vertexing cuts. The parameters u mg and /i a ' slg (signal events) are included as 
free parameters in the global fit, and are assumed to be the same for all B° peaking background 
sources. For B + peaking background components, the z^s and ^s are fixed to the values extracted 
from a previous unbinned maximum-likelihood fit to the tagged and untagged At distributions of 
B + data used as control samples, described in Sec. 4. We neglected and ^x a ^ for combinatorial 
background sources. 
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7.3 Signal and background characterization 



The function Tree (C) i n Eq. (30) describes the signal or background probability of observing a 
particular value of £■ It satisfies 

/Cm ax 
dC£«(C) = l, (37) 
7>min j 

where [Cmim Cmax] is the range of wies/AE values used for analysis. 

For -Bfl av and B CPK o events, the t«es shape is described with a single Gaussian distribution 
for the signal and an ARGUS parameterization for the background [18]. Based on these fits, an 
event-by-event signal probability prec g (wiEs) can be calculated for each tagging category a and 
sample rec. As we do not expect signal probability differences between B° and B°, the ttt-es fits 
were performed to B^ av and -Bfj av events together. Due to the lack of statistics and the high purity 
of the sample, the to-es fits to the B° — ► ip(2S)Kg and B° — ► XciK® samples were performed 
without splitting by tagging category. The component fractions TrecirriEs) are then given by 



^ee ig (^Es) = [l " / r a e'r k J PrfcH^Es) 
•T^rVEs) = f^PrfcH^Es) 

^cVes) = [l-^ s c ig (m ES )]/ r a e f (38) 
where (3 indexes the combinatorial, prompt and non-prompt, background components, 

E/- =1 " (39) 



The fraction /" e ' pcak of the signal Gaussian distribution is due to backgrounds that peak in the 
same region as the signal, determined from Monte Carlo simulation [18]. The estimated contri- 
butions are (1.5 ± 0.6)%, (0.28 ± 0.11)%, (1.8 ± 0.6)%, (1 ± 3)% and (3.5 ± 1.4)% for the 5 flav , 
J/^K° S {K° s -► 7T+7T-), J/0#° {K° s -► 7r°7r°), ^{2S)K° S and X ciK° s channels, respectively. A com- 
mon peaking background fraction is assumed for all tagging categories within each decay mode. 
We take a common prompt fraction for all tagging categories for each B CPK o decay channel in- 
dependently. Due to the higher statistics of the Bq &v sample and the significant differences in 
the background levels for each tagging category, /^o Prompt = /^o Pr ° mP is allowed to vary with the 

tagging category. Note that the parameters of the JF" e 'c lg (mEs) functions, determined from the set 
of separate unbinned maximum-likelihood fits to the itt-es distributions, are kept fixed in the global 
fit. 

For B CPK o events the background level is much higher, with significant non-combinatorial 
components, therefore requiring special treatment [18]. A binned likelihood fit to the AE spectrum 
in the data is used to determine the relative amounts of signal and background from B — >■ J/ip X 
events and events from a misreconstructed J/ip — > l + £~ candidate (non-J/0 background). In these 
fits, the signal and B — > J/ipX background distributions are obtained from inclusive J/tjj Monte 
Carlo, while the non- J/ip distribution is obtained from the J/0 dilepton mass sideband. The Monte 
Carlo simulation is also used to evaluate the channels that contribute to the B — > J/ipX background. 
Due to differences in purity and background composition, the fit is performed separately for each 
reconstruction type (EMC and IFR) and lepton type (J/0 — ► e + e~ and J/0 — > The different 
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inclusive J/ij; backgrounds from Monte Carlo are then renormalized to the J/ij; background fraction 
extracted from the data. The fractions are adjusted for lepton-tagged and non-lepton-tagged events 
in order to compensate for the observed differences in flavor tagging efficiencies in the J/ip sideband 
events relative to the -Bfl av and inclusive J/ip Monte Carlo. In addition, some of the decay modes 
in the inclusive J/ip background have particular CP content. The PDF can then be formulated as 



+ 



j=JH> * 



^CPK° 



(3 



(40) 



As the J/ip lepton type is not expected to influence the AE shape, the PDFs are used without 
regard to lepton type. The AE PDFs are used separately for EMC and IFR K® type, and they 
are grouped for J/ipK® (signal), J/ipK® background, J/ipX background (excluding J/ip Kg) and 
non- J/ip . 



7.4 Signal and background structure 

For signal events, a common set of mistag, At resolution function, and v and fi a parameters for 
all samples are assumed. This assumption is supported by extensive Monte Carlo studies. Peaking 
backgrounds originating in B° decays are assumed to have the same parameters as the signal. For 
B + peaking backgrounds we still assume the same resolution function as for signal, but the mistag 
and v, fj, a parameters are fixed to the values extracted from the previous unbinned maximum- 
likelihood fit to the E> + data. For combinatorial background components we use an empirical 
description of the mistags and At resolution, allowing various intrinsic time dependencies. In the 
nominal fit we assume prompt and non-prompt components, the non-prompt component being a 
pure exponential dependence. As discussed in Sec. 10, fits with oscillatory, CP/T, CPT/CP and 
DCKM structure have been also performed to evaluate possible systematic biases. A common set 
of mistags and At resolution parameters, independent of the signal, is assumed for the non-J/fy 
background in the B CPK o sample and for the prompt and non-prompt background components in 
the i?fl av and B CPK o samples. In this case, the parameters Aw a and w°[ are fixed to zero, and 
the resolution model uses the core and outlier Gaussian distributions. The fraction of prompt and 
non-prompt component and the exponential constant of the non-prompt component in the non- J/ip 
background are fixed to the values obtained from an external fit to the time distribution of the 
J/ip dilepton mass sideband. The nominal fit neglects i/i and /i aj for combinatorial background 
sources. 



7.5 Free parameters for the nominal fit 

The ultimate aim of the fit is to obtain simultaneously sgn(ReAcp)Ar/r, \q/p\, (Re\cp/\^cp\) ^ ez 
and Imz, assuming Rqp = 1. The parameters lm\cp/\\cp\ and Am are also floated to account for 
possible correlations and to provide an additional cross-check of the measurements, comparing our 
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values with the BABAR sin 2/3 result based on the same data sample [14] and recent Ato measure- 
ments. The average B° lifetime, tb, is fixed to the PDG value, 1.542 ps [19]. As a cross-check we 
also performed fits with Rcp and tb allowed to vary. All these physics parameters are by construc- 
tion common to all samples and tagging categories, although the statistical power for determining 
each parameter comes from a particular combination of samples or At dependences, as discussed 
in Sec. 2. 

The doubly-CKM-suppressed parameters are necessarily small and difficult to determine. In 
particular, ReAfl av occurs only multiplied by other small parameters. As a result, we have chosen to 
fix ReAfl av = 0. We fit for the parameter ImAfl av /|Afl av |, but fix |^4fl av /Aiav| = 0.02. We regard these 
two variables as formal entities, without the connection implied by taking ReAfl av = 0. We vary 
separately ImAfl av /|Afl av |, but keep |Afl av | = |Afl av | \p/q\ 2 - Thus there are two free parameters, plus 
one fixed magnitude. We treat the tagging DCKM similarly. We assign a systematic error by vary- 
ing the magnitudes by 100% and scanning all possible combinations of the phases (Sec. 10). This 
strategy provides, for the current data sample size, the optimal trade-off for all measured parame- 
ters between statistical and systematic uncertainties originating in the ambiguities of doubly-CKM 
suppression in the B tag , and to a lesser extent the -Bfl av . 

The total number of parameters that are free in the fit is 58, of which 36 parameterize the 
signal: physics parameters (4), cross-check physics parameters (2), single effective imaginary parts 
of the doubly-CKM-suppressed phases (4), resolution function (10), mistags (11) and differences in 
the fraction of B° and B° mesons that are tagged and reconstructed (5). The At distributions, the 
asymmetries and the physics parameters sgn(ReAcp)Ar/T, \q/p\, (ReAcp/|Acp|) Rez, Imz and the 
cross-check parameter ImAcp/|Acp| were kept hidden until the analysis was finished. However, the 
parameter Am, the residual At distributions and asymmetries, the statistical errors and changes 
in the physics parameters due to changes in the analysis were not hidden. 

8 Analysis results 

We extract the parameters sgn(ReAcp)Ar/T, \q/p\, (ReAcp/|Acp|) Rez and Imz, ImAcp/|Acp|, 
Am, the parameters for doubly-CKM-suppressed decays, the signal mistag, resolution function and 
v, fi a parameters and the empirical background parameters with the likelihood function described 
in Sec. 7. In Table 3 we list the signal yields in each tagging category after vertexing requirements. 
The purities (estimated in the region toes > 5.27 GeV/c 2 for non-B CPK o samples and \AE\ < 10 
MeV for B CPK o events), averaged over tagging categories, are 82%, 94% and 55%, for -Bfl av , B CPK o 
and B CPK ^ candidates, respectively. The fitted signal mistag and resolution function parameters 
are shown, respectively, in Tables 4 and 5. The values of the asymmetries in reconstruction and 
tagging efficiencies are summarized in Table 6. There is good agreement with the asymmetries 
extracted with the counting-based approach outlined in Appendix A. Note that in the counting 
technique the time-dependence of UnTagged events is not used, therefore decreasing significantly 
the statistical power of the measurement of AT/T. 

The values of sgn(ReAcp)Ar/r, \q/p\, (ReAcp/|Acp|) Rez and lmz extracted from the fits are 
given in Table 7. The fitted effective doubly-CKM-suppressed decay parameters are also indi- 
cated. All these results can be compared to those obtained when the fit was repeated assuming 
CPT invariance. The significant change of the effective doubly-CKM-suppressed decay parameters 
between the two fits is due to the large correlation of these parameters with the CPT violating 
parameter lmz. The fitted value of Ato agrees with recent S-Factory measurements [1, 2, 3, 4], 
and remains unchanged between the two fits. The fit result for ImAcp/|Acp| when we assume CPT 
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Table 3: Signal event yields, obtained from the ttjes fits f° r the -Bflav and B CPK o samples and 
multiplying by the signal fraction in the \AE\ < 10 MeV interval for the B CPK o sample, after 
vertexing requirements. 





-Bflav 


B CPK% 


B CPK° Ti 


Tag 




W 


Tot 




w 


Tot 


B* 


w 


Tot 


Lepton 


1478 


1419 


2897 


96 


98 


194 


35 


35 


70 


Kaon I 


2665 


2672 


5337 


154 


175 


329 


74 


65 


139 


Kaon II 


3183 


2976 


6159 


181 


188 


369 


85 


66 


151 


Inclusive 


3197 


3014 


6211 


184 


172 


356 


78 


72 


150 


UnTagged 




10423 






585 






260 





Table 4: Signal tagging efficiencies and mistag parameters for each tagging category a as extracted 
from the nominal maximum-likelihood fit. Uncertainties are statistical only. 



Tagging category 


t a (%) 


< 


^slopc 


Aw a 


Lepton 


9.4 ±0.2 


0.026 ± 0.007 


(fixed) 


-0.012 ±0.012 


Kaon I 


17.2 ±0.3 


0.020 ± 0.020 


0.13 ±0.04 


-0.027 ±0.013 


Kaon II 


19.9 ±0.3 


0.159 ±0.024 


0.07 ±0.04 


-0.042 ±0.013 


Inclusive 


19.9 ±0.3 


0.265 ±0.025 


0.07 ±0.04 


-0.029 ±0.013 


UnTagged 


33.6 ±0.6 









Table 5: Signal At resolution function parameters as extracted from the nominal maximum- 
likelihood fit. Uncertainties are statistical only. 



Parameter 


Fitted value 


Parameter 


Fitted value 


Score 


1.25 ±0.04 


•S'tail 


5.7 ±0.8 


rLepton 

Ocore 


0.02 ±0.07 


<5tail 


-1.5 ±0.5 


rKaonl 

"core 


-0.27 ±0.05 


/tail 


0.034 ±0.010 


rKaonll 

"core 


-0.32 ± 0.04 


Cout 


8 ps (fixed) 


^Inclusive 

"core 


-0.30 ± 0.04 


"out 


ps (fixed) 


rUnTagged 

"core 


-0.28 ±0.03 


/out 


0.0003 ±0.0012 



invariance agrees with our sin 2(5 measurement based on the same data set [14]. When we allow for 
CPT violation, IvaXcp /\\cp\ increases by ±0.012 with unchanged statistical errors. The statistical 
correlation coefficients among all physics and cross-check physics parameters are shown in Table 
8. Table 9 shows the top five statistical correlations of the physics parameters with any other free 
parameter in the global fit. 

Figures 2 and 3 show the At distributions of the signal candidates (to-es > 5.27 GeV/c 2 for £>fl av 
and B CPK o and \AE\ < 10 MeV for B CPK o samples). The points correspond to data. The curves 
correspond to the projections of the nominal likelihood fit weighted by the appropriate relative 
amounts of signal and background. The background contribution is indicated with the shaded 
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Table 6: Values of the signal B°B° differences in reconstruction (y) and tagging (n a ) efficiencies as 
extracted from the nominal maximum-likelihood fit. The results are compared with those obtained 



Parameter 


Nominal fit 


Counting-based method 


V 


O.Olf ±0.008 


0.007 ± 0.008 


..Lepton 

H- 


0.024 ± 0.022 


0.029 ± 0.042 


..KaonI 

H' 


-0.022 ±0.017 


-0.022 ± 0.029 


..KaonI I 

H 1 


0.014 ±0.016 


0.004 ± 0.027 


Inclusive 

H- 


0.014 ±0.016 


0.025 ± 0.027 



Table 7: Physics parameters results, from the global nominal fit and when we assume CPT invari- 
ance. The free single effective doubly-CKM-suppressed decay parameters are also indicated. Errors 
are statistical only. 



Parameter 


Nominal fit results 


Fit results assuming CPT invariance 


sgn(ReA C p)Ar/r 


-0.008 ± 0.037 


-0.009 ± 0.037 


\Q/P\ 


1.029 ±0.013 


1.029 ±0.013 


(Re\ C p/\Xcp\)Rez 


0.014 ± 0.035 




Imz 


0.038 ± 0.029 




ImAtag/ Atag 


1.5 ± 1.2 


0.5 ± 1.0 


ImAtag/ A tag 


-0.1 ± 1.2 


0.8 ±1.0 


ImAflav/ Aflav 


2.3 ±1.1 


1.4 ±0.9 


ImAflav/ A flav 


-0.6 ±1.1 


0.1 ±0.9 



Table 8: Correlation among all the physics parameters as extracted from the simultaneous 
maximum-likelihood fit to the -Bfl av and Bcp samples. 





sgn(ReA C p)Ar/r 


\Q/P\ 


ImAcp/ \Xcp\ 


(ReXcp/\X C p\)Rez 


Imz 


Am 


-1.3% 


-2.8% 


-5.6% 


7.0% 


-0.2% 


sgn(ReA C p)Ar/r 




11.0% 


0.4% 


-7.9% 


-1.8% 


\q/p\ 






-1.0% 


-2.4% 


-1.1% 


ImAcp/ \Xcp\ 








-10.9% 


17.4% 


(ReX C p/\Xcp\)Rez 










-3.4% 



area. 



9 Cross-checks and validation studies 

We have used data and Monte Carlo samples to perform validation studies of the analysis 
technique. Monte Carlo tests include studies with parameterized fast Monte Carlo as well as full 
GEANT4 [17] simulation samples. Checks with data were performed with control samples where 
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Table 9: Top five correlations of the physics parameters with any other free parameter of the global 
maximum-likelihood fit. 



Phvmrs 'n^raTTiptPT' 


Pa ra rn pf pr 


(lorrplatiori ( %) 


Am 


c 

"-Hail 


— 10.4 




1 \ ^CP \ ) J^ez 

...KaonI 

w 


o.o 
-5.8 




...Kaonll 
^slope 
cUnTagged 

Ucorc 


/I V 

—4. f 
4.2 


sgn(ReA C p)Ar/r 




11.0 




{Ke\cp/\Acp\) tiez 


— ( .y 




V 


0.0 




A „.,KaonI 

law 


— 3.8 




A „,,KaonI 


Q /I 

— o.4 


W/P\ 


V 

A „,, Kaonll 


65.1 

99 ^ 
— Zz.O 




A „,,KaonI 

LIU/ 


99 4 




A Tnrl n^i vp 
Za7/? 


—15.5 




..Kaonll 


3.2 


ImX C p/\X C p\ 


Imz 


17.4 




T™ \ /I \ 1 

lmA tag/ | A tag| 


14.4 




ImAflav/IAflavl 


—6.2 




ImAtag/|Atag| 

KaonTT 
^Kaonll 


-5.5 
4.5 


(ReAcp/|Acp|)Rez 


lm\ C p/\Xcp\ 


-10.9 




Imz 


q /i 




Lepton 


2.2 




ImAtag/ Atag| 
Kaonll 


-2.0 
-1.6 


Imz 


ImAtag/ A tag 


61.6 




ImAtag/|Atag| 
Au ,KaonI 


-56.6 
-8.0 




Au ,KaonII 


-6.1 




Au ,Inclusive 


-3.5 



no Ar, CP jT or CPT/CP effects are expected. Other checks were made by analyzing the actual 
data sample, but using alternative tagging, vertexing and fitting configurations. 

9.1 Monte Carlo simulation studies 

An extensive test of the fitting procedure was performed with fast parameterized Monte Carlo 
simulations, generating 300 experiments with a sample size and composition corresponding to that 
of the data. The mistag rates and At distributions were generated according to the model used in 
the likelihood function. The nominal fit was then performed on each of these experiments. Each 
experiment used the set of m,Es{A-E) and a At values observed in the non-K^ (K®) sample. The 



27 



Oh 


400 


in 




o 


300 


t/3 




+- > 

rt 




> 


200 








100 





in 

olOOO 

t/3 
+- > 

W 500 







: a) Mixed events _| 


BaBar 




\ B° tags 




+ B° tags 




^^^■**«— 


b) Unmixed events ±* 


\ B° tags 




+ B° tags 




■ — -'- -'- -i- -i- - 



10 







10 
At (ps) 



Figure 2: The At distributions for (a) -Bfl av mixed and (b) -Bfl av unmixed events with a B° tag and 
with a B° tag in the signal region, toes > 5.27 GeV/c 2 . The solid (dashed) curves represent the 
fit projection in At based on the individual signal probabilities and event-by-event At uncertainty 
for B° (B°) tags. The shaded area shows the background contribution to the distributions. 

rms spread of the residual distributions for all physics parameters (where the residual is defined as 
the difference between the fitted and generated values), was found to be consistent, within 10%, 
with the mean (Gaussian) statistical errors reported by the fits. Moreover, from these experiments 
it was verified that the asymmetric 68% and 90% confidence intervals provide the correct coverage. 
The mean values of the residual distributions were in all cases consistent with no measurement bias. 
A systematic error due to the precision of this study was assigned to each physics parameter. The 
statistical errors on all the physics parameters (Table 7) and the calculated correlation coefficients 
among them (Tables 8 and 9) extracted from the nominal fit to the data were consistent with the 
range of values obtained from these experiments: We found that 24% of the fits result in a value 
of the negative log-likelihood that is less (better) than that found in data. 

In addition, large samples of signal and background Monte Carlo events generated with a full 
detector simulation were used to validate the measurement. The largest samples were generated 
with zero values of AT/T, \q/p\ — 1 and z, but additional samples were also produced with rela- 
tively large values of these parameters. Other values (including those measured in the data) were 
generated with reweighting techniques. The signal Monte Carlo samples were split into samples 
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Figure 3: The At distributions for (a) B CPK o and (b) B CPK o events with a B° tag and with a 
B° tag in the signal region, uies > 5.27 GeV/c 2 for B CPK o candidates and \AE\ < 10 MeV for 
B CPK o^ events. The solid (dashed) curves represent the fit projection in At based on the individual 

signal probabilities and event-by-event At uncertainty for B° (B°) tags. The shaded area shows 
the background contribution to the distributions. 



comparable to the actual data set, keeping the relative sizes of signal -Bfj av , B CPK o and B CPK o 
samples as observed. To check whether the selection criteria, or the analysis and fitting procedure, 
introduced any bias in the measurements, the nominal fit (signal only) was then applied to these 
experiments. The small combinatorial background in these signal samples was rejected, using only 
events in the signal region. Fits to the pure signal physics model, using the true At distribution 
and true tagging information, were also performed. The means of the residual distributions from 
all these experiments for all the physics parameters are consistent with zero, confirming that there 
is no measurement bias. The rms spreads are consistent with the average reported errors and with 
the estimated errors in the fit to data. A systematic error is assigned to each physics parameter due 
to the limited Monte Carlo statistics from this test. The effect of backgrounds has been evaluated 
by adding an appropriate fraction of background events to the signal Monte Carlo sample and 
performing the fit. The Bqp background samples were obtained either from simulated B — > J/ipX 
events or AE sidebands in data, while the B^v backgrounds were obtained from generic BB Monte 
Carlo. We find no evidence of bias in any of the physics parameters. 
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9.2 Cross-checks with data 

We fit subsamples selected by choosing a tagging category and running period. Fits using only 
the B° -► L>M-X+ or B° -► J/ip K*°(K + ir-) specific channels for £ flav , and 5 CPX o or B CPK o 
only for i?cp> were also performed. We found no statistically significant differences in the results 
for the different subsets. We also varied the maximum values of |At| and a At accepted between 5 
and 30 ps, and 0.6 and 2.2 ps, respectively. Again, we did not find statistically significant changes 
in the fitted values of the physics parameters. 

In order to verify that the results are stable under variation of the vertexing algorithm that 
is used in the measurement of At, we used alternative (less powerful) methods [18]. In order 
to reduce statistical fluctuations due to different events being selected, the comparison between 
the alternative and nominal methods was performed using only the common events. Observed 
variations were consistent with anticipated statistical fluctuations and with the systematic error 
assigned to the resolution function (see Sec. 10). 

The stability of the results under variation of the tagging algorithm was studied by repeating 
the fit using the tagging algorithm used in Ref. [18]. The algorithm used in that analysis had an 
effective tagging efficiency, Q = Y^ a t a (1 ~~ 2w a ) 2 , about 7% lower than the one used here. The 
variation observed in the physics parameters is consistent with the statistical differences. 

The average B° lifetime was fixed in the nominal fit to the PDG value [19]. This value was 
obtained by averaging measurements based on flavor eigenstate samples and assuming negligible 
effects from AT/T, \q/p\ and CPT violation. Measurements that do not use tagged events are 
not affected by \q/p\ and CPT violation, but are by a non-zero value Ar/r at second order, as 
discussed in Sec. 2. Therefore we do not expect significant effects from the fixed average B° lifetime. 
However, to check the consistency of the result, the fit was repeated with tb left free. The resulting 
tb was about two standard deviations below the nominal value assumed in our analysis, before bias 
corrections and taking into account the statistical error from the fit and the present tb uncertainty. 
The changes of the physics parameters were within statistical errors. Nevertheless, a systematic 
error was assigned using the variation of each physics parameter repeating the fit with tb fixed 
after a change of two standard deviations of the present uncertainty (±0.032 ps). 

Similarly, fits with Rqp floated were performed. The resulting Rcp value was consistent with 1 
(the fixed nominal value), within one standard deviation (statistical only). The changes observed 
in the physics parameters were consistent with their statistical uncertainties. Systematic errors due 
to fixing Rcp at unity were set by changing Rcp by twice the statistical uncertainty determined by 
floating it (±10%). The resulting variation in each parameter was taken as the systematic error. 

The robustness of the fit was also tested by modifying the nominal PDF normalization, as 
described by Eq. (33), so that the analysis was insensitive to the relative amount of B° and B° 
tagged events. As a consequence, the statistical power of \q/p\ was dramatically reduced, since 
the sensitivity of this parameter comes largely from the differences in time-integrated B° and B° 
rates. In addition, the fit was also performed assuming an independent set of resolution function 
parameters for each tagging category. In all cases the results are compatible within statistical 
differences with the nominal fit results. Finally, the tagging efficiencies i aj were alternatively 
determined from the -Bfl av sample and assumed to be the same for all samples (rather than to use 
the estimate from each sample separately, as in the nominal fit). The change in the values of the 
physics parameters was found to be negligible. 

Control samples in data from B + decays (treated in a fashion analogous to that described in 
Sec. 4) have also been used to validate the analysis technique, since in these samples we expect zero 
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values for AT/T, \q/p\ and z. For the -Bfl av sample we used the B + — > D 7r + decay channel, and 
for the Bcp sample the charmonium B + decays. Due to the absence of mixing and CP violation 
in these samples, the check was performed fixing Am = and \q/p\ = 1 in the -Bfl av sample, and 
Am = 0.489 ps _1 and ImXcp /\Xcp\ = in the Bqp sample, fitting only for sgn(ReAcp)Ar/r, 
(ReAcp/|Ac*p|) Rez and Imz. No statistically significant deviations from zero were observed. 

10 Systematic uncertainties 

We estimate systematic uncertainties with studies performed on both data and Monte Carlo 
simulation samples. A summary of the non-negligible sources and results is shown in Table 10. In 
the following, the individual contributions are referenced by the lettered lines in this table. 

10.1 Likelihood fit procedure 

Several sources of systematics due to the likelihood fit procedure are considered. We first 
include the results from the tests performed using the fast parameterized Monte Carlo (a) and the 
full GEANT4 signal Monte Carlo events (b), as described in Sec. 9.1. We take the larger of the 
observed bias (mean of the residual distributions) and its statistical error as the systematic error. No 
corrections are applied to the central values extracted from the fit to the data since no statistically 
significant bias is observed. Note that the GEANT4 simulation addresses the underlying assumption 
that the fit properly accounts for residual differences in the mistag, resolution function, v, fi a 
parameters for -Bfl av , B CPK o and B CPK ^ samples and differences in At resolution for correct and 
wrong tags [18]. 

We also consider the impact on the measured physics parameters of the asymptotic PDF nor- 
malization. The effect is evaluated by varying the fitted values using a normalization in the range 
defined by the At cut. Finally, the fixed tagging efficiencies are varied within their statistical 
uncertainties. The two contributions are found to be negligible. 

10.2 At resolution function 

The resolution model used in the analysis, consisting of the sum of three Gaussian distributions, 
is expected to be flexible enough to accommodate the actual resolution function. To assign a 
systematic error to the assumption of this model, we use the alternative model described in Sec. 6, 
with a Gaussian distribution plus the same Gaussian convolved with one exponential function, for 
both signal and background. The results for all physics parameters obtained from the two resolution 
models are consistent within statistical uncertainties. However, we assign the difference of central 
values as a systematic uncertainty (c). 

In addition a number of parameters that are integral to the determination of At are varied 
according to known uncertainties. The PEP-II boost, estimated from the beam energies, has an 
uncertainty of 0.1% [16]. The absolute z-scale uncertainty has been evaluated to be less than 0.4%. 
This estimate was obtained by measuring the beam pipe dimensions with scattered protons and 
comparing to optical survey data. Therefore, the boost and z-scale systematics are evaluated by 
varying by ±0.6% the reconstructed At and a/s.t (d). The uncertainty on the beam spot, which 
is much wider than it is high, is taken into account by moving its vertical position (the direction 
most valuable in vertexing) by 20 and 40 fim and increasing the vertical dimension by 30 and 60 
fim (e) . Finally, the systematic uncertainty due to possible S VT internal misalignment is evaluated 
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Table 10: Summary of systematic uncertainties on the measurement of sgn(ReAcp)Ar/r, \q/p\, 
(ReAcp/|Acp|) Rez and Imz. 



Systematics source 


sgn(ReA C p)Ar/r 


\q/p\ 


(ReAcp/|A C p|)Rez 


Imz 


Likelihood fit proced 


ure 






(a) Parameterized MC test 


0.003 


0.001 


0.003 


0.003 


(b) GEANT4 MC test 


0.005 


0.007 


0.004 


0.016 


At resolution function 


(c) Res. funct. parameterization 


0.007 


0.001 


0.008 


0.003 


(d) z scale and boost 


0.003 


0.001 


0.002 


< 0.001 


(e) Beam spot 


0.008 


0.002 


0.001 


0.011 


(f) SVT alignment 


0.006 


0.001 


0.001 


0.011 


(g) Outliers 


0.002 


< 0.001 


< 0.001 


< 0.001 


Signal properties 


(h) Average B° lifetime 


0.004 


0.001 


0.004 


< 0.001 


(i) Direct CP violation 


0.002 


0.004 


0.001 


0.003 


(j) Doubly-CKM-suppressed decays 


0.008 


0.004 


0.032 


0.006 


(k) Residual charge asymmetries 


0.005 


0.006 


0.004 


0.006 


Background properties and structure 


(1) Signal probability 


0.002 


0.001 


0.002 


0.001 


(m) Fraction of peaking background 


< 0.001 


< 0.001 


0.004 


< 0.001 


(n) At structure 


0.002 


0.001 


0.001 


0.001 


(o) Ar/CPT/CP/T/Mixing/DCKM 


0.001 


0.002 


0.002 


< 0.001 


(p) Residual charge asymmetry 


< 0.001 


0.001 


< 0.001 


< 0.001 


(q) K\ specific systematics 


0.004 


< 0.001 


0.004 


0.003 


Total systematics 


0.018 


0.011 


0.034 


0.025 



by applying a number of possible misalignment scenarios to a sample of simulated events and 
comparing the values of the fitted physics parameters from these samples to the case of perfect 
alignment (f). 

Fixing the width and bias of the outlier component, respectively to 8.0 and 0.0 ps, introduces 
systematic errors. To estimate the uncertainty we add in quadrature the variation observed in the 
physics parameters when the bias changes by ±5 ps, the width varies between 6 and 12 ps and the 
outlier distribution is assumed to be flat (g). 
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10.3 Signal properties 

As described in Sec. 9.2, the uncertainty from fixing the average B° lifetime has been evaluated 
by moving its central value by ±0.032 ps (h), twice the current uncertainty [19]. Possible direct 
CP violation in the Bqp sample was taken into account by varying Rcp by ±10% (i). 

Systematics from doubly-CKM-suppressed decays arise due to uncertainties in ReAtag, ReA ta g, 
ReAfl av and ReAfl av - In order to evaluate this contribution, samples of parameterized Monte Carlo 
tuned to the data sample, with all possible values of the doubly-CKM-suppressed phases were 
generated, assuming a single hadronic decay channel contributing to the -Btag and to the -Bflav 
The generation was made for maximal values of |A ta g/A;ag| an d l^flav/^flavl) assuming a 100% 
uncertainty on its estimate based on the elements of the CKM matrix, ~ 0.02 [19]. For the 
Lepton tagging category, largely dominated by semileptonic B decays, we assumed |A tag | to be 
zero. Using the fit results from all these samples, we evaluated the larger of the offset with respect 
to the generated value and its statistical uncertainty, for each possible configuration of phases. The 
systematic error assigned was the largest value among all configurations (j). This is the dominant 
source of systematic uncertainty for the measurement of (ReAcp/|Acp|) Rez, primarily due to the 
DCKM effects in the tagging B meson. Similar studies performed with more than one hadronic 
channel indicated the destructive interference among the different decay modes, proving that our 
prescription to assign the systematics assuming a single effective channel is conservative. 

Charge asymmetries induced by a difference in the detector response for positive and nega- 
tive tracks are included in the PDF and extracted together with the other parameters from the 
time-dependent analysis. Thus, they do not contribute to the systematic error, but rather are 
incorporated into the statistical error at a level determined by the size of the -Bfl av data sample. 
Nevertheless, in order to account for any possible and residual effect, we assigned a systematic 
uncertainty as follows. We reran the B reconstruction, vertexing and tagging code after killing ran- 
domly and uniformly (no momentum or angular dependence) 5% of positive and negative tracks 
in the full Monte Carlo sample. This 5% is on average more than a factor three larger than the 
precision with which the parameters v and ^ a have been measured in the data. The half difference 
between the results obtained for positive and negative tracks is assigned as a systematic error (k). 

10.4 Background properties and structure 

The event-by-event signal probability prtc g (w-es) f° r -Bflav and B CPK o samples was fixed to the 
values obtained from the tbes fits. We compared the results from the nominal fits to the values ob- 
tained by changing one sigma up and down all the tties distribution parameters, taking into account 
their correlations. This was performed simultaneously for all tagging categories, and independently 
for the i?fiav and B CPK o samples. Alternatively, we also used a flat signal probability distribution: 
events belonging to the sideband region (mES<5.27 GeV/c 2 ) are assigned a signal probability of 
zero, while we gave a signal probability equal to the purity of the corresponding sample to signal 
region events (m.Es>5.27 GeV/c 2 ). The differences among fitted physical parameters with respect 
to the default method were found to be consistent within the statistical differences. We determined 
the systematic error due to this parameterization by varying the signal probability by the statis- 
tical error in the purity. The final systematic error was taken to be the larger of the one-sigma 
variations found for the two methods (1). The uncertainty on the fraction of peaking background 
was estimated by varying the fractions according to its uncertainty separately for the -Bfl av sample 
and each B CPK o decay mode (m). The effective r\cp of the B CPK o peaking background, assumed 
to be zero in the nominal fit, was also varied between ±1 and —1 and found to be negligible. 
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Another source of systematic uncertainty originates from the assumption that the At structure 
of the combinatorial background in the m-ES sideband region is a good description of the structure 
in the signal region. However, the background composition changes slightly as a function of m^s, 
since the fraction due to continuum production slowly decreases towards the B mass. To study this 
effect, we first varied the lower edge of the t«es distributions from 5.20 GeV/c 2 to 5.27 GeV/c 2 , 
simultaneously for the Ba av and B CPK o samples, observing a good stability of the result. We also 
split the sideband region in seven equal slices each 10 MeV/c 2 wide and used each of these ranges 
to perform a standard fit. The quadratic sum of the extrapolation to the B mass region and the 
error on it was assigned as systematic uncertainty (n). 

As described in Sec. 7, the nominal likelihood fit assumes that there is no Ar, CP/T, CPT/CP, 
mixing and doubly-CKM-suppressed decays content in the combinatorial background components 
(-Bflav and B CPK o samples) and in the non-J/ifj background (B GPK o sample). To evaluate the effect 
of this assumption we repeated the fit but now assuming non-zero values of Ar, \q/p\, z, ImAcp and 
Am, varying rjcp of the background by ±1. The check was performed by introducing in the PDF 
an independent set of physics parameters and assuming maximal mixing and CP violation (Am 
and lm\cp/\\cp\ fixed to 0.489 ps _1 [19] and 0.75 [14], respectively). Doubly-CKM-suppressed 
decay effects were included assuming the maximal values of | A ta _ g / A ta _ g \, |^4fl a v/^4flav| and scanning 
all the possible values of the B° and B° phases for B^ av and .Btag- The systematic uncertainty was 
evaluated simultaneously for all of these sources (o) . 

The uncertainty due to the B + lifetime has been evaluated by moving the central value according 
to the current uncertainty [19]. It was found to be negligible. The B + mistags and the differences in 
the fraction of B + and B~ mesons that are tagged and reconstructed were varied according to their 
statistical errors as obtained from the fit to the B + data. They were found also to be negligible. 
Uncertainties from charge asymmetries in combinatorial background components (neglected in the 
nominal fit) were evaluated by repeating the fit with a new set of v and /i Q parameters. The 
measured values of v and \x a are found to be compatible with zero and the variation of the physical 
parameters with respect to the nominal fit is assigned as systematic error (p). 

For the B GPK o channel [18], the signal and non-J/i/; background fractions are varied according 
to their statistical uncertainties as obtained from the fit to the AE distribution. We also vary 
background parameters, including the J/ipX branching fractions, the assumed r]cp, the AE shape 
and the fraction and effective lifetime of the prompt and non-prompt non-J/^/; components. The 
differences observed between data and Monte Carlo simulation for the K® angular resolution and 
for the fractions of B° — ► J/ipK® events reconstructed in the EMC and IFR are used to evaluate 
a systematic uncertainty due to the simulation of the K® reconstruction. Finally, an additional 
contribution is assigned to the correction applied to Lepton events due to the observed differences 
in flavor tagging efficiencies in the J/ip sideband relative to -Bfl av and inclusive J/tp Monte Carlo. 
Conservatively, this error was evaluated comparing the fit results with and without the correction. 
The total B CPK o specific systematics is evaluated by taking the quadratic sum of the individual 
contributions (q). 

10.5 Summary of systematic uncertainties 

All individual systematic contributions described above and summarized in Table 10 are added 
in quadrature. The dominant source of systematic error in the measurement of (ReAcp/| Acp|) Rez 
is due to our limited knowledge of the doubly-CKM-suppressed decays, which also contributes sig- 
nificantly to the other measurements. The limited Monte Carlo statistics are a dominant source 
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of systematics for \q/p\, Imz and to a lesser extent to sgn(ReAcp)Ar/T. Residual charge asym- 
metries also dominate the systematics on \q/p\. Our limited knowledge of the beam spot and 
SVT alignment also reflects significantly on Imz and sgn(ReAcp)Ar/T. The systematic error on 
sgn(ReAcp)Ar/r also receives a non- negligible contribution from our understanding of the resolu- 
tion function. The systematic uncertainties on sgn(ReAc*p)Ar/T and \q/p\ when CPT is assumed 
to be a good symmetry were evaluated similarly, and found to be, respectively, ±0.020 and ±0.012. 

11 Summary and discussion of results 

The conventional analysis of mixing and CP violation in the neutral B meson system neglects 
possible contributions from several sources that are expected to be small. These include the dif- 
ference of the lifetimes of the two neutral B meson mass eigenstates, the CP- and T-violating 
quantity \q/p\ — 1, which is proportional to Im(Ti2/Mi2), and potential CPT violation. To measure 
or extract limits on these quantities requires the full expressions for time dependence in mixing and 
CP violation and consideration of systematic issues that might mimic fundamental asymmetries 
we seek to measure, like detector charge asymmetries, different resolution function for positive and 
negative At, and doubly-CKM-suppressed decays from both fully reconstructed final flavor states 
and non-leptonic tagging states. 

Our analysis of approximately 31,000 fully reconstructed flavor eigenstates and 2600 CP eigen- 
states sets new limits on the difference of decay widths of B° mesons and on the CP, T, and CPT 
violation intrinsic to B°B° mixing. The six independent parameters governing mixing (Am, AT/T), 
CPT /CP violation (Rez, Imz) and CP /T violation (ImAcp, \q/p\) are extracted from a single fit of 
both fully reconstructed CP and non-CP events, tagged and untagged. This provides the sensitiv- 
ity required to separate the small effects we seek from asymmetries in detector response and from 
potentially obscuring correlations in the decays of the two B mesons. The preliminary results are 

sgn(ReA C p)Ar/r = -0.008±0.037(stat.) ± 0.018(syst.) [-0.084,0.068] , 

\q/p\ = 1.029±0.013(stat.) ±0.011(syst.) [ 1.001,1.057], 

(ReAcp/|A C p|)Rez = 0.014±0.035(stat.) ± 0.034(syst.) [-0.072,0.101], 

Imz = 0.038±0.029(stat.) ± 0.025(syst.) [-0.028,0.104] . 

The values in square brackets indicate the 90% confidence-level intervals. When estimating the 
limits we also evaluated multiplicative contributions to the systematic error, adding in quadrature 
with the additive systematic uncertainties. Assuming CPT invariance the results are 

sgn(ReAcp)Ar/r = -0.009±0.037(stat.) ± 0.020(syst.) [-0.087,0.069] , 
\q/p\ = 1.029±0.013(stat.) ± 0.012(syst.) [ 1.000,1.058]. 

The parameters Am and ImAcp/|Acp| are allowed to float, so that recent U-Factory Am results 
[1, 2, 3, 4] and our sin 2(3 analysis based on the same data sample [14] provide a cross-check. The 
value of the CP/T- violating parameter ImAcp/|Acp| increases by +0.012 when CPT violation is 
allowed in the fit. 

The results are consistent with Standard Model expectations and with CPT invariance. To 
date, these are the best limits on the difference of decay widths of B° mesons and the strongest 
test of CPT invariance outside the neutral kaon system [12]. The limit on CP and T violation in 
mixing is independent of and consistent with our previous measurement based on the analysis of 
inclusive dilepton events [20]. Fig. 4 shows the results in the (|(//p| — 1, \z\) plane, comparing them 
to the BABAR measurement of \q/p\ made with dileptons and to the Standard Model expectations. 
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All the other results are also consistent with previous analyses [4, 6, 7, 13]. While the Standard 
Model predictions for Ar and \q/p\ — 1 are still well below our current limits and no CPT violation 
is anticipated, higher precision measurements may still bring surprises. 




Figure 4: Favored regions at 68% confidence level in the (\q/p\ — 1, \z\) plane determined by this 
analysis ("Hadronic") and by the BABAR measurement of the dilepton asymmetry [20]. The axis 
labels reflect the requirements that both CP and T be violated if \q/p\ / 1 and that both CP 
and CPT be violated if \z\ ^ 0. The region shown for this analysis is constrained to lie within 
the physical region \z\ > and is chosen to exclude the maximum range of alternative hypotheses 
given the a priori low probability of obtaining \z\ =0 due to phase-space considerations. The 
dilepton measurement constrains \q/p\ without assumptions on the value of \z\. The Standard 
Model expectation is obtained from Ref. [21]. 
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A Efficiency asymmetries 

The use of untagged data is essential to determining the asymmetries in the tagging and re- 
construction efficiencies. To indicate how the various samples enter we provide a simple example 
using only time-integrated quantities. In practice we use a time-dependent analysis, which gives 
better precision because it uses more information. Suppressing the tag category a, the signal or 
background component, j, and writing the reconstruction efficiencies as r = r J B0 , r = r^ and the 

tagging efficiencies as t = t°^a , t = , Eq. (36) reads 

r — r 



r + r 



" - TTV <41> 

Using the numbers of signal events that are tagged and have a reconstructed B° (X), those tagged 
and having a B° (Y), those untagged with a reconstructed B° (Z) and finally those untagged 
with a reconstructed B° (W) we can determine the required asymmetries [18]. To see this, note 
that if the total number of B°B° pairs is N, and neglecting Ar, \q/p\ and z corrections, there are 
N u = N(l + [1/(1 + a$)])/2 unmixed events (i.e. B°B°) and N m = N(l - [1/(1 + x 2 d )})/2 mixed 
events (i.e. B°B° or B°B°), where x,i = rgAm, so 

X = rtN m /2 + riN u /2 
Y = rtN m /2 + rtN u /2 
Z = r{l-t)N m /2 + r(l-t)N u /2 

W = r(l-t)N m /2 + r(l-t)N u /2 . (42) 
Setting U = X + Z and V = Y + W, we find 

,- u ~ v u - (l | x2) (y/v)-(x/u) 

U ~U + V ' ^- {1 + Xd) (Y/V) + (X/U) ■ (43) 

Corrections to these equations have to be applied due to non-zero values of Ar, \q/p\ — 1 and z. 
The use of untagged events is essential to the determination of v and jjl. 
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